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We first discuss how the longstanding confusion in the literature concerning one-loop 
quantum corrections to 1 + 1 dimensional soli tons has finally been resolved. Then we 
use’t Hooft and Veltman’s dimensional regularization to compute the kink mass, and 
find that chiral domain wall fermions, induced by fermionic zero modes, lead to spon¬ 
taneous parity violation and an anomalous contribution to the central charge such that 
the BPS bound becomes saturated. On the other hand, Siegel’s dimensional reduction 
shifts this anomaly to the counter terms in the renormalized current multiplet. The 7 • j 
superconformal anomaly is located in an evanescent counter term, and imposing su¬ 
persymmetry, this counter term induces the same anomalous contribution to the central 
charge. Next we discuss a new regularization scheme: local mode regularization. The 
local energy density computed in this scheme satisfies the BPS equality (it is equal to 
the local central charge density). In an appendix we give a very detailed account of the 
DHN method to compute soliton masses applied to the supersymmetric kink. 


1 Introduction 

Quantum corrections to solitons were of great interest in the 1970’s and 1980’s 
PUB 01. and again in the last few years, due to the present activity in quantum field 
theories with dualities between extended objects and pointlike objects. Dashen, 
Hasslacher, and Neveu [JTp, in a 1974 article that has become a classic, computed 
the one-loop corrections to the mass of the bosonic kink in A field theory and to 
the bosonic soliton in sine-Gordon theory. For the latter, there exist exact analytical 
methods associated with the complete integrability of the system, authenticating the 
perturbative calculation. Our work here uses general principles but focuses on the 
kink, for which exact results are not available. Dashen et al. put the object (clas¬ 
sical background field corresponding to kink or to sine-Gordon soliton) in a box 
of length L to discretize the continuous spectrum, and used mode number regular¬ 
ization (equal numbers of modes in the topological and trivial sectors, including 
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the zero mode in this counting) for the ultraviolet divergences. They imposed peri¬ 
odic boundary conditions (PBC) on the meson field which describes the fluctuations 
around the trivial or topological vacuum solutions, and added a logarithmically di¬ 
vergent mass counter term whose finite part was fixed by requiring absence of tad¬ 
poles in the trivial background. They found for the one-loop correction to the kink 
mass 


= X) \hu, n - £ i&4°> + SM = -Urn < 0 , (1) 

where m is the mass of the meson in the trivial background and SM the counterterm 
induced by renormalizing m. This result remains unchallenged. 

The supersymmetric (susy) case, as well as the general case including fermions, 
proved more difficult. The action reads 

£ = -\ - \ u " - ’ (2) 

where — |f/ 2 = —1(<^ 2 — /i 2 /A) 2 , the meson mass is m = /i\/2, and c = 1 for 
supersymmetry. Dashen et al. did not publish the fermionic corrections to the 
soliton mass, stating “The actual computation of [the contribution to] M W [due to 
fermions] can be carried out along the lines of the Appendix. As the result is rather 
complicated and not particularly illuminating we will not give it here” (page 4137 

of©). 

Several authors have since performed the calculation of M (l > for the susy kink, 
and found different answers. It became clear that the answers depended on the 
choice of boundary conditions (BC) for the fluctuation fields, more precisely on the 
BC for the fermions. Moreover, it also became clear that one obtained different 
answers if one used different regularization schemes. At present these issues are 
believed to be fully understood as follows. 

Boundary conditions: Boundary conditions distort fields near the boundary. 
This distortion creates spurious boundary energy which should be subtracted from 
the total energy in order to obtain the true mass of the kink. There are several ways 
to avoid the spurious boundary energy 

(i) one may first compute the energy density £ (x) and then integrate over a re¬ 
gion which contains the kink but stays away from the boundaries 0; 

(ii) one may average overs sets of BC such that in the average the boundary en¬ 
ergy cancels [0. One such set of BC for fermions which has been studied 
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in detail consists of periodic BC (ip ± (—L/2) = antiperiodic BC 

(^ ± (—L/2) = —'0 ± (L/2)), twisted periodic BC (^(—L/ 2) = i> T (L/ 2)), 
and twisted antiperiodic BC 2) = — -0 =F (L/2))Q 

(iii) one may choose a set of BC which have no boundary. By this cryptic state¬ 
ment we mean BC which put the system on a circle (more precisely a Mobius 
strip) such that the system becomes translationally invariant and one cannot 
identify a point where the boundary is present [|]]. In principle such BC 
could still lead to delocalized (homogeneously spread out) boundary energy, 
but this does not occur [^]. By using the Z 2 symmetry ip K (—x) = —<Pk(%) of 
the ki nk background, one such set of BC has been identified to be the twisted 
(anti)periodic BC in the kink sector; 

(iv) one may first consider a background which contains both a kink (K) and an 
antikink (K) with periodic BC, and then divide the answer for the mass of this 
compound KK system by 2 [||, |5|]. (Putting a kink next to an antikink, there 
is a small cusp in the background where the kink is joined to the antikink, 
but for large distances the effect of this cusp can be neglected. One can also 
find an exact solution which is everywhere smooth and has periodic BC (a 
“sphaleron”) but this involves transcendental functions.) In fact, if one begins 
with periodic BC for the fermions in this KK system, one finds that the mode 
solutions have either twisted periodic or twisted antiperiodic BC in between. 

Regularization schemes: Several well-known regularization schemes have been 
applied to the calculation of the quantum kink mass and the quantum central charge. 
To regulate the various sums over zero-point energies one has used: mode number 
cutoff, energy-momentum cutoff, heat-kernel techniques, (^-function techniques, ’t 
Hooft-Veltman’s dimensional regularization, Siegel’s dimensional regularization 
(“dimensional reduction”). To regulate Feynman graphs, one has used higher- 
space-derivative regularization with factors (1 — d * 2 /M 2 ) (this regularization of the 
kinetic terms but not the interactions preserves susy, although it breaks Lorentz in¬ 
variance^ and again dimensional regularization. It has turned out that the reason 
some of these schemes give incorrect answers is that they were applied incorrectly: 
one naively applied the rules which had been developed for a trivial background 
to the kink background. After proper modification, these schemes all now yield 
the same answers. It is of some interest (and useful for avoiding errors in future 
calculations) to point out the required modifications of all these schemes. In the 

'Strictly speaking, these BC should be called even and odd rather than periodic and antiperiodic. 

2 Because the anticommutator of two supersymmetry charges never produces a Lorentz generator, 
it is possible to preserve supersymmetry while breaking Lorentz symmetry. 


3 



following, however, we concentrate on discussing in detail the two variants of di¬ 
mensional regularization as well as a newly proposed method to study the local 
energy distribution of the quantum mass, local mode regularization 00]. 

The ’t Hooft-Veltman dimensional regularization can be employed in a susy 
preserving manner by embedding the minimally susy kink in 1 < d < 2 spatial 
dimensions. This leads to new physics, namely spontaneous parity breaking and 
chiral domain wall fermions, which provide a new explanation [llj for the origin 
of the anomalous contribution [|]] to the central charge of the susy kink. Siegel’s 
dimension reduction, on the other hand, where d < 1, obtains this anomaly from 
an evanescent counter term to the superconformal current, which gives rise to an 
anomalous nonconservation at the quantum level of the conformal version of the 
central-charge current [If]. 


2 Dimensional regularization and reduction 


2.1 One-loop bosonic kink mass 


Probably the most elegant regularization scheme to avoid the difficulties of mode 
regularization in a finite box and the possibility of boundary energy is dimensional 
regularization by embedding the 1 + 1 dimensional ki nk in n = d + 1 dimensions as 
a domain wall. 


As has been shown in Ref. [|12|], this reproduces correctly the one-loop quantum 
mass of the bosonic 1+1 dimensional kink, as well as the surface tension of the 
higher-dimensional kink domain walls [ |T3| ] . 

By analytic continuation of the number of extra transverse dimensions (d — 1) 
of a kink domain wall, no further regularization is needed. Denoting the momenta 
pertaining to the extra dimensions by t and reserving k for the momentum along 
the kink, i.e. perpendicular to the kink domain wall, the energy of the latter per 
transverse volume L d ~ 1 is obtained from summing/integrating zero-point energies 
according to 


mm 

L d ~ x 


m 


3A 

1 



d d ~H 

( 2 ? r ) d — 1 



00 rlkd d ~ 1 f 

j- V k 2 + l 2 + m 2 S' K (k) + 5M 

oo 


( 3 ) 


where the discrete sum is over the normalizable states B of the 1 + 1-dimensional 
kink with energy lvb, and the integral is over the continuum part of the spectrum. 
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The spectrum of fluctuations for the 1+1-dimensional kink is known exactly 
[f]~4|] . It consists of a zero-mode, a bound state with energy uj 2 R /rri 2 = 3/4, and 
scattering states in a reflectionless potential for which the phase shift 5 K (k) = 
—2 arctan(3 mk/ ( m 2 — 2k 2 )) in the ki nk background provides the difference in the 
spectral density, 6' K {k), between kink and trivial vacuum. 

In a “minimal” renormalization scheme where tadpoles cancel but Z\ = 1, one 
has 

SM = ^ d 2 j dk{k 2 + m 2 )^+ (4) 

r(- 5 )(47T) 2 Jo 

yielding (with x = k/m ) 


ilfW 

L d - 1 


m 3 T(^)m d r 1 /3\ 2 

3-^ r(-i)(4vr)V 12 V 4 7 

f dx(x 2 + l)^ 2 


-1 

4x 2 + 1 


-(ri- 1) 


( 5 ) 


Here the first term within the braces is the contribution from the bound state with 
nonzero energy, and the second is the result of combining the last two terms in ( 3]). 
In the limit d —> 1, which corresponds to the 1+1 dimensional kink, one obtains 


„ ail m3 ( m 3 m\ 

d = 1= 3x + U7^~^) , 


( 6 ) 


reproducing the well-known DHN result |JT[]. It is interesting to note that it is the 
last term in (|S[) that would be missed in a sharp-cutoff calculation (see Ref. O) 
and that it now arises from the last term in the square brackets of (|5|). 

Eq. (5j) is also valid for d —* 2 where it gives the surface tension of a 2+1 
dimensional kink domain wall; for higher dimensions one has to include also a 
renormalization of the coupling A. All these results are in agreement with those 
obtained previously by other methods [|T5|]. 


2.2 One-loop susy kink mass 

Dimensional regularization is more delicate in susy theories. To preserve susy, 
one should normally consider Siegel’s dimensional regularization by dimensional 
reduction [|T6], [l7j]. However, it is also possible to preserve susy by embedding the 
susy kink in dimensions <2 + 1. 

Embedding the susy kink in 2+1 dimensions gives a domain wall centered about 
a one-dimensional string on which the fermion mass vanishes (since U'{^k) °c 
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p K vanishes at the center of the kink). The total energy M of the domain wall is 
infinite but the energy density M/L is finite; as a result there is strictly speaking 
no zero mode in 2+1 dimensions associated with translational invariance. Indeed, 
the zero mode of the ki nk is only normalizable in 1+1 dimensions, but one can 
construct eigenfunctions in 2+1 dimensions which are products of zero modes in 
1 + 1 dimensions and plane waves in the orthogonal direction(s) (along the domain 
wall). 

The 2+1 dimensional case is different also with respect to the discrete symme¬ 
tries of (|2j). In 2+1 dimensions, y 5 = y^y 2 = ±1 corresponding to the two 
inequivalent choices available for y 2 = ±r 1 (in odd space-time dimensions the 
Clifford algebra has two inequivalent irreducible representations). Therefore, the 
sign of the fermion mass (Yukawa) term can no longer be reversed by 0 —> y 5 0 
and there is no longer the Z 2 symmetry p —> —p, 0 —> y 5 0. 

What the 2+1 dimensional model does break spontaneously is instead parity, 
which corresponds to changing the sign of one of the spatial coordinates. The La- 
grangian is invariant under x m —> —x m for a given spatial index m — 1,2 together 
with p —> —p (which thus is a pseudoscalar) and 0 —> y m 0. Each of the trivial 
vacua breaks these invariances spontaneously, whereas a kink background in the re¬ 
direction with <Pk(~‘ x 1 ) = —c Pk (0 1 ) is symmetric with respect to x l -reflections, 
but breaks x 2 = y reflection invariance. 

This is to be contrasted with the 1 + 1 dimensional case, where parity (x l —> 
—x 1 ) can be represented either by 0 —> y°0 and a true scalar ip —> p or by 0 ► y 3 0 
and a pseudoscalar ip —» —p. The former leaves the trivial vacuum invariant, and 
the latter the ground state of the kink sector. 

In what follows we shall consider the quantum corrections to both, the mass 
of the susy kink and the tension of the domain string, together. We again use a 
minimal renormalization scheme, where inclusion of the fermionic tadpole loop 
simply replaces the prefactor 3 in (|j) by (3 — 2). 

In a Majorana representation of the Dirac matrices in terms of the usual Pauli 
matrices r k with y° = — ir 2 , y 1 = r 3 , y 2 = r 1 (added for d = 2), and C = t 2 so 
that -0 = (Jy) with real 0 + (x, t ) and 0 _ (x, t), the equations for the bosonic and 
fermionic normal modes with frequency u ; and longitudinal momentum £ (nonzero 
only when d = 2) in the kink background p = px read 


[~d 2 x + u ' 2 + uu"\n = (w 2 

(d x + f/')0 + + i{u + f)0 - 

( d x — U')^~ + i(tv — £)0 + 


t 2 )v, 

0, 

0 . 


(V) 

( 8 ) 

( 9 ) 


Acting with ( d x — U') on (S) and eliminating 0 as well as p' = —V shows that 


6 


0 + satisfies the same equation as the bosonic fluctuation rj. Compared to ip + , the 
component 0~ has a continuous spectrum whose modes differ by an additional 
phase shift 6 = — 2 arctan(m/ k) when traversing the kink from x\ = — oc to x i = 
■foo, which is determined only by U'(<pk(xi = ±oo)) = ±m. Correspondingly, 
the difference of the spectral densities of the 0 + -fluctuations in the kink and in the 
trivial vacuum equals that of the ^-fluctuations, whereas that of 0“-fluctuations is 
obtained by replacing 5' K —>► 5' K + O'. 

In the sum over zero-point energies for the one-loop quantum mass of the kink 
(when d = 1), 


M — M cL + - iv B ~ ^ — 2 UJp ~ a; ^) ^ 

one thus finds that the bosonic contributions from the continuous spectrum are can¬ 
celed by the fermionic contributions]] except for the additional contribution involv¬ 
ing Q'(k) in the spectral density of the -0“ modes. 

The discrete bound states cancel exactly, apart from the subtlety that the fermionic 
zero mode should be counted as half a fermionic mode [5|]. In strictly 1+1 dimen¬ 
sions, the zero modes do not contribute simply because they carry zero energy, and 
for d > 1, where they become massless modes, they do not contribute in dimen¬ 
sional regularization. 

In a cutoff regularization in d = 2, as we shall further discuss below, they in fact 
do play a role. Remarkably, the half-counting of the fermionic zero mode for d = 1 
has an analog for d — 2 where the bosonic and fermionic zero modes of the kink 
correspond to massless modes with energy |az| = \l\. From (|8]) and (g) one finds 
that the fermionic ki nk zero mode 0+ oc cp' K , -0“ = 0 is a solution only for uj = -\-£. 

It therefore cancels only half of the contributions from the bosonic kink zero mode 
which for d = 2 have u = ±L For d — 2 one thus finds that the fermionic zero 
mode of the kink corresponds to a chiral (Majorana-Weyl) fermion on the domain 
wall (string) in 2+1 dimensions. O0I 

In dimensional regularization, however, the kink zero modes and their massless 
counterparts for d > 1 can be dropped, and the energy density of the susy domain 
wall reads 


M l 1 ) m 3 1 
TF 1 = 3A ~ 4 


dkd d ~H 

(2vr) d 


\/k 2 + £ 2 + m 2 9\k) + SM, 


(ID 


J This cancellation could be however incomplete for certain boundary conditions in global mode 
regularization. 

4 Choosing a different sign for 71 reverses the allowed sign of £ for these fermionic modes and 
thus their chirality (with respect to the domain string world sheet). This corresponds to the other, 
inequivalent representation of the Clifford algebra in 2+1 dimensions. 
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where 


( 12 ) 


0\k) = 


2m 


k 2 + 


m “ 


With 5M = \5M the logarithmic divergence in the integral in ( |TT) ) as d — > 1 
gets cancelled. A naive cut-off regularization at d = 1 would actually lead to a total 
cancellation of the /('-integral with the counter term 5M, giving a vanishing quantum 
correction in renormalization schemes with A = Ao- In dimensional regularization 
there is now however a mismatch for d ^ 1 and a finite remainder in the limit d—>\ 
proportional to (d — l)r (—(d — l)/2). The final result reads [ [131 ] 


L d ~ l 


m 

3A 


m d 2 
( 47 t )^ d 


r( 


3 — d 
2 ^ 


(13) 


In view of the discussion of the central charge below, it is instructive to write 
the above finite remainder that dimensional regularization leaves behind for d —* 1 
in the form 


MW m 3 _ 1 fdkd d ~Hf 

IT 11 ~ 3A “ ~4 J (2vr ) d 7° ' > 


(14) 


which is obtained by combining the integral in ([TT]) with the integral representation 
of the counter term (1/3 of the r.h.s. of (Q)). Evidently, the nonvanishing result is 
entirely due to the momenta in the extra d — 1 dimensions of a kink domain wall. 

In the literature, at least to our knowledge, only the case of a supersymmetric 
kink (d — 1) has been considered and dimensional regularization reproduces the 
result obtained before by Refs. [j<|. [20|. |6. [TT). j4j, [22[ |. 

However, a (larger) number of papers have missed the contribution —m/ (27r), 
mostly because of the (implicit) use of an inconsistent energy-cutoff scheme [ j23[ , 
000 or have obtained different answers because of the use of boundary con¬ 
ditions that accumulate a finite amount of energy at the boundaries [ [27[ [T5[ ] . The 
former result is however now generally accepted and, in the case of the super-sine- 
Gordon model (where the same issues arise with the same results) in agreement 
with S-matrix factorization P28j]. 

A new result, which follows from ( fT3| ) and which will play a role for the discus¬ 
sion of central charges in the next section, is the nonvanishing one-loop correction 


M { dl 2 _ 

L 87t 


(15) 


for the surface tension of the minimally susy kink domain wall in 2+1 dimensions. 

In Ref. [129] the correct susy kink mass has also been obtained by employing a 
smooth energy (momentum) cutoff, the necessity of which becomes apparent, as in 
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the purely bosonic case, by considering the 2+1 dimensional domain wall. Using 
a naive cutoff for d — 2 one finds quadratic divergences which cancel only upon 
inclusion of the zero modes (which become massless modes in 2+1 dimensions). As 
we have discussed above, unlike the other bound states, these do not cancel because 
the fermionic zero mode becomes a chiral fermion on the domain-string world- 
sheet and thus cancels only half of the bosonic zero (massless) mode contribution, 
yielding 



Vk 2 + P + m 2 —: -, 

k z + m 


A^oo r M(£_ 
~Jo 2tt12 


l t 1 

— arctan — 1 
7r mJ 


'0 


00 dim 
7r 27T 


y/k 2 + P + m 2 


(16) 


which is however still linearly divergent. Smoothing out the cutoff in the k- integral 
does pick an additional (and for d — 1 the only) contribution —rnj (27t), which is 
now necessary to have a finite result for d = 2. This finite result then reads 


M, 


(i) 


d =2 


1 f°° dt ( m 

— — m — i arctan — 

7T In 27T V £ 


m 

87T 


(17) 


in agreement with the result obtained above in dimensional regularization. 


2.3 Susy algebra and its quantum corrections 

2.3.1 Dimensional regularization 


The susy algebra for the 1+1 and the 2+1 dimensional cases can both be covered by 
starting from 2+1 dimensions, the 1+1 dimensional case following from reduction 
by one spatial dimension. 

In 2+1 dimensions one obtains classically [|3U[| 


{Q a ,Q p } = 27(7 M )7m, (M = 0,1,2) 

= 2i( 1 °H + 1 1 (P x + Z y )+ 1 2 (P y ~Z x )) a p , 

where we separated off two surface terms Z m in defining 

P m = [ d d xV m , Vm = fidmV - 


7 — / 7 

I u/ i 


Zm = U{(p)d m (f = d m W{ip) 


(18) 


(19) 

( 20 ) 
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with W(p) = f d^p U(ip). 

Having a kink profile in the ^-direction, which satisfies the Bogomolnyi equa¬ 
tion d x pK = —U(pk), one finds that with our choice of Dirac matrices 


Q ± = j d 2 x[(p =f dyip)^ + (d x p ±U(p))^}, (21) 

{Q ± ,Q ± } = 2(H±(Z X -P y )), (22) 

and the charge Q + leaves the topological (domain-wall) vacuum p = p K , -0 = 0 
invariant. This corresponds to classical BPS saturation, since with P x = 0 and 
P v = 0 one has |Q + , Q + | = 2 (H + Z r ) and, indeed, with a kink domain wall 
Z x /L d ~ 1 = W(+v) - W(-v) = -M/L d ~\ 

At the quantum level, hermiticity of Q ± implies 

{s\H\s) > |(s|P,|s)| = |<s|(P s -Z*)|s)|. (23) 


This inequality is saturated when 


Q + |s) = 0 


(24) 


so that BPS states correspond to massless states PmP m = 0 with P y = M for a 
kink domain wall in the x-dircction pT|], however with infinite momentum and en¬ 
ergy unless the ^/-direction is compact with finite length L. An antikink domain wall 


has instead Q~ 
broken. 


s') = 0. In both cases, half of the supersymmetry is spontaneously 


Classically, the susy algebra in 1+1 dimensions is obtained from ( |T8l ) simply by 
dropping P y as well as Z y so that P x = P x . The term 7 2 Z X remains, however, with 
7 2 being the nontrivial 7 s of 1+1 dimensions. The susy algebra simplifies to 


{Q ± : Q ± } = 2(H±Z), {Q + ,Q~} = 2P X (25) 


and one has the inequality 

(s\H\s) > |(s|Z|s)| (26) 

for any quantum state s. BPS saturated states have Q + |s) = 0 or Q~ |s) = 0, corre¬ 
sponding to kink and antikink, respectively, and break half of the supersymmetry. 

In a kink (domain wall) background with only nontrivial x dependence, the 
central charge density Z x receives nontrivial contributions. Expanding Z x around 
the ki nk background gives 

= Ud x p K - =d x p K + d x (Ur)) + ^djU'r) 2 ) + 0(r] 3 ). (27) 
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where only the part quadratic in the fluctuations contributes to the integrated quan¬ 
tity at one-loop ordeiQ. However, this matches precisely the counter term 5M from 
requiring vanishing tadpoles. Straightforward application of the rules of dimen¬ 
sional regularization thus leads to a null result for the net one-loop correction to 
(Z x ) in the same way as found in Refs. [ j24j p5j [13116|] in other schemes. 

On the other hand, by considering the less singular combination (H + Z x ) and 
showing that it vanishes exactly, it was concluded in Ref. [j2l|l that (Z x ) has to com¬ 
pensate any nontrivial result for (H ), which in Ref. [|2Tj| was obtained by subtracting 
successive Bom approximations for scattering phase shifts. In fact, Ref. [ pT[ ] explic¬ 
itly demonstrates how to rewrite (Z x ) into —(H), apparently without the need for 
the anomalous terms in the quantum central charge operator postulated in Ref. [Q]. 

The resolution of this discrepancy is that Ref. [ ]2Ij ] did not regularize (Z x ) and 
therefore the manipulations needed to rewrite it as — (H) (which eventually is regu¬ 
larized and renormalized) are ill-defined. Using dimensional regularization naively 
one in fact obtains a nonzero result for (H + Z x ), apparently in violation of susy. 

However, dimensional regularization by embedding the kink as a domain wall 
in (up to) one higher dimension, which preserves susy, instead leads to 

(H + Z x — P y ) = 0, (28) 


i.e. the saturation of (fZ3j), as we shall now verify. 
The bosonic contribution to ( P y ) involves 


1 

2 


(rjdyV + d y r]rj) 



(29) 


where the 0k(x) are the mode functions of the fluctuation field operator //. The £- 
integral factorizes and gives zero both because it is a scale-less integral and because 
the integrand is odd in L 

The fermions on the other hand turn out to give nontrivial contributions: The 
mode expansion for the fermionic field operator reads 


p = 


V’o + 
4>o = 


d 


(2tt) 

' d ^ 1 


bl, (c.c.) 


H yf dk 


(30) 


where ipo is the fermionic zero-mode lifted to a Majorana-Weyl domain wall fermion, 

5 But this does not hold for the central charge density locally [J|, M|. 
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and s k = (d x + U')^. This leads to 

(Py) = 


4 


dk l 
2n 2 lo 


2 XT y 

1 r d d ~ l 

2 J ( 2 tt ^- 

1 /' d d ~ l f 

-2Swr' te( - t)W+ 

1 f d d ~ x t v-£' dk 

+ 2 ] W) 


[(a; + £)\4>k\ 2 + (a; — f)|sfc| 2 ] 




2t \2 (W +|Sil )+ 2^ (W 


- |s»| 2 ) • (31) 


From the last sum-integral we have separated off the contribution of the zero mode 
of the kink (the chiral domain wall fermion for d > 1). The contribution of the latter 
no longer vanishes by symmetry, but the ^-integral is still scale-less and therefore 
put to zero in dimensional regularization. The first sum-integral on the right-hand 
side is again zero by both symmetry and scalelessness, but the final term is not: 
The (-integration no longer factorizes because u = \/k 2 + i 2 + m 2 , and leads to 
a nonvanishing result, which, as one can show 0L is identical to the finite net 
contribution in (H). For the integrated quantities, this equality can be seen by 
comparing with ( [74] ) upon using that f dx(\4> k \ 2 — |s fc | 2 ) = 9'{k). 

So for all d < 2 we have BPS saturation, (H) = \ (Z x — P y ) |, which in the limit 
d —> 1, the susy kink, is made possible by a nonvanishing (P y ). The anomaly in 
the central charge is seen to arise from a parity-violating contribution in d = 1 + e 
dimensions which is the price to be paid for preserving (minimal) supersymmetry 
when going up in dimensions to embed the susy kink as a domain wall. 

To summarize, in 2+1 dimensions, we have P y = P y — Z x and (P y )[ = ( H ), 
where P and Z were defined in (|l8j). Classically, this BPS saturation is guaranteed 
by Z x alone. At the quantum level, however, the quantum corrections to the latter 
are cancelled completely by the counter term from renormalizing tadpoles to zero. 
All nontrivial corrections come from the “genuine” momentum operator P y , and 
are due to having a spontaneous breaking of parity. 

In the limit of 1+1 dimensions, because 7 2 |d= 2 +i = 7 5 |d=i+i> one has to make 
the identification Z = Z x — P y . For Z x , one again does not obtain net quantum cor¬ 
rections. However, the expectation value (P y ) does not vanish in the limit d —> 1, 
although there is no longer an extra dimension. The spontaneous parity violation in 
the 2+1 dimensional theory, which had to be considered in order to preserve susy, 
leaves a finite imprint upon dimensional reduction to 1+1 dimensions by provid¬ 
ing an anomalous additional contribution to (Z x ) balancing the nontrivial quantum 
correction ( H). 
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2.3.2 Dimensional reduction 


We now show how the central charge anomaly can be recovered from Siegel’s ver¬ 
sion of dimensional regularization [ fT5| , [T7| ] where n is smaller than the dimension 
of spacetime and where one keeps the number of field components fixed, but lowers 
the number of coordinates and momenta from 2 to n < 2. At the one-loop level 
one encounters 2-dimensional 5" coming from Dirac matrices, and n-dimcnsional 
6” from loop momenta. An important concept which is going to play a role is that 

of the evanescent counterterms [|32|] involving the factor jd", where 8” = 8" — 8” 
has only e = 2 — n nonvanishing components. 

Consider now the supercurrent j lt = — (0ip + U (ip) YpPP. In the trivial vacuum, 
expanding into quantum fields yields 

, 7 /i = - + U'(v ) rj + ^U"(v) r 7 ^ + -^=^ 2 7 ^, (32) 


where v = /i/\/X Only matrix elements with one external fermion are divergent. 
The term involving U"(v)r ) 2 in ( [32] ) gives rise to a divergent scalar tadpole that is 
cancelled completely by the counter term ()jr (which itself is due to an 77 and a ip 
loop). The only other divergent diagram is due to the term involving I 77 in ( [32[ ) and 
has the form of a U-sclfcncrgy. Its singular part reads 


(0|j M b) dlv = iU"(v) / dx 
Jo 


d n K 


tint 


(27r) n [Y 2 + p 2 x( 1 — x) + m 2 } 


-u(p). 


(33) 


Using 8 V = 8" — 5" we find that under the integral 

- -^)7a = -« 2 7/x - -« 2 ^7 a 
^ n ^ n n ^ 


so that 


<<w“ v = 


d.» _ U"(v) si 


7 a u(p). 


(34) 


27r e 

Hence, the regularized one-loop contribution to the susy current contains the evanes¬ 
cent operator 


•div = U"iv) 

J LL 


7a'0- 


(35) 


^ 2t r e 

This is by itself a conserved quantity, because all fields depend only on the n- 
dimensional coordinates. The renormali z ed susy current f™' = Jp, — j'f 1; is thus still 
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conserved^] but from the evanescent counter term it receives a nonvanishing contri¬ 
bution to 7 • f™' which appears in the divergence of the renormalized conformal- 
susy current jfi"' 

S'W’-W = -Tjf’ = —!>■ (36) 

(There are also nonvanishing nonanomalous contributions to d p ( fij p ) because our 
model is not conformal-susy invariant at the classical level.) 

Ordinary susy on the other hand is unbroken; there is no anomaly in the diver¬ 
gence of f™' ■ A susy variation of j p involves the energy-momentum tensor and the 
topological central-charge current ( /; . according to 

= - 2 T/ 7.6 - 2C M 7 5 0 (37) 

where classically = e /u/ Ud l 'ip. 

At the quantum level, the counter-term = —j p lv induces an additional con¬ 
tribution to the central charge current 

C“ = \~ e ^ dPU> (38) 

which despite appearances is & finite quantity: using that total antisymmetrization 
of the three lower indices has to vanish in two dimensions gives 

(39) 

and together with the fact the U' only depends on n-dimensional coordinates this 
finally yields 

C“ = d " U ' < 4 °) 

in agreement with the anomaly in the central charge as obtained previously.|] 

We emphasize that C /t itself does not require the subtraction of an evanescent 
counterterm. The latter only appears in the susy current j /1 , which gives rise to a 
conformal-susy anomaly in yfcj^. A susy variation of the latter shows that it forms a 
conformal current multiplet involving besides the dilatation current T /u /i M and the 
Lorentz current T^ u x p e up also a current 

jfj) 1 = x p e pv C- (41) 

6 Note also that ( 35|) does not change the susy charges Q = f dxjo if one assumes that S p has 
only spatial components. Furthermore, recall that conserved currents do not renormalize. 

7 It would be interesting to study further the infrared/ultraviolet connection for this anomaly. 
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We identify this with the conformal central-charge current, which is to be distin¬ 
guished from the ordinary central-charge current 

The anomalous contribution to the ordinary central charge is thus understood 
as the additional nonconservation at the quantum level of the conformal central- 
charge current (pTT|). (Additional, because the model is not conformally invariant so 
that there is already nonconservation at the classical level.) This finally answers the 
question: what kind of anomaly corresponds to the anomalous contribution to the 
central charge? 

3 Local Casimir energy for solitons 

We have seen in the introduction that there is a problem with the regularization of 
the zero-point energies by means of mode number cutoff (equal numbers of modes 
in each sector, with careful counting of zero modes): it includes spurious boundary 
energy. On the other hand, the principle of mode regularization seems natural, so 
the question arises whether we can devise a mode number cutoff scheme without the 
unwanted boundary energy. This almost automatically leads to a new regularization 
scheme for Casimir sums, called local mode regularization. Given that each mode 
determines a mode function 4> n (x, t ) (or f> n (x, t ) for the fermions) normalized such 
that for a large box with volume L the t ) become at large \x\ a plane wave 
with unit strength (the corrections to these plane waves are of order L -1 / 2 ), we can 
introduce a concept of local mode density p{x) in the ki nk sector (and p <iv> (x) in the 
trivial sector) as follows: p(x) = f* n {x)f n {x). To regulate such sums we would 
like to again cut off the sum over n at a large number N. 

The kink mass contains the difference of the energy sum \cj) n (x)\ 2 ^u n in the 
kink sector and the energy sum \ fn\x)\ 2 \un^ in the trivial sector. The problem 
is thus how to relate the regularization in one sector to that in the other sector. The 
most straightforward method would be to include the same number of terms in each 
sum, just as in the case of global mode number regularization. However in this way 
the sums would only indirectly take the presence of the kink into account (through 
the inequality of the to n and uJn } ). 

We now formulate a principle which we have not yet been able to prove that it is 
equivalent to other principles, or that it preserves supersymmetry, but which gives 
correct answers for the kink mass and supersymmetric kink mass, and which is so 
simple that it deserves further study. Namely we require that the regulated mode 
densities in both sectors are equal. The function Yln=o is a function of 

N or equivalently of A = 2nN/L, but for large L we can interpolate it to become a 
function of a continuous variable A. Similarly, p® becomes a continuous function 
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of A. Since L~ l J \<p k \ 2 dx = 1 counts each mode once, it may seem natural to 
also use \(j> k {x)\ 2 to count modes locally. However note that (<p(x)<p(x)) contains 
22 \4> k (x)\ 2 while the energy density contains 22 2 \^k- The choice to 

use 22 \(f> k (x)\ 2 to define a regularization is perhaps natural, certainly more natural 
than for example 22 | 4 , but we have not proven that 22 \<Pk{x) 1 2 is the correct 

object. 

If the density p(x) is cut off at a large A, the density p (0) (x) should be cut off 
at a A + A A = 2n A® such that p A (x) is equal to Pa+aa( x )- F ar awa Y from the 
kink all modes are plane waves, so for large \x\ one expects A A to vanish, but near 
the kink AA will be nonvanishing. This implies that AA is x-dependent, and the 
principle of local mode regularization takes the following form 

PA Or) =PA+ AA( x)( a; )- ( 42 ) 

The regulated energy densities in the kink and the trivial sector, given by e = 
22 \4>n{x)\ 2 \u n and e® = 22 \(j)n\x)\ 2 \un\ will then i n general be different if 
the regulated densities p(x) and p (0 ) (x) are equal. 

It is now straightforward to calculate the local Casimir mass of a soliton. It is 
given by 

ecasOr) = e(x) - e (0) (x) = 

1 f A fjlf 1 / , A+AA(x) J u 1 

-u B <t>%(x) + 2j ^m,x)\ 2 -u-2^ --u: + 5M{x). (43) 

The bound state yields a zero-point energy 4 cu B and has (normalizable) mode func¬ 
tion 4>b{x), while the continuous spectrum in the kink sector consists of plane waves 

x) = exp ikx. We rewrite this expression such that it is manifestly finite 

£Cas(%) = \u B <t> 2 B {x) 

+ {2 £ A (|0(x, k)\ 2 - 1) + SM(x)\ - 4A!A. (44) 

The last term is the ’’anomaly”, it appears here as a term of the form A/A because 
AA(x) is proportional to 1/A as well see presently. 

In the kink sector (f>(k, x) can be given explicitly. From the explicit form one 
finds that it can be expressed in terms of the wave functions of the discrete spectrum 
as follows 

2,/m 2 — u 2 

|0(A;,x)| 2 - 1 =-^0 2 (x)-^ - —. (45) 

j 3 
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(For the kink j refers to the zero mode with Uj — 0 and the bound state with 
Uj = 3/4 m 2 .) This formula seems to be new and we interpret it as a local version of 
the completeness relation. Integration over k yields the usual completeness relation 

/ °° (Ik 

— {\(j)(k,x)\ 2 - 1 } + 4%{x) + <t> 2 B {x) = 0, (46) 

but the local version allows us to evaluate AA as 


AA(x) 


dk (l — | <t>(k, tc)| 2 ) 


3 rrr 


4A cosh r -f- 


0{ A 


—2\ 



(47) 


The local counter term SM(x) is of course equal to the term proportional to d/r 
in the energy density.0 We can now substitute all these relations and find then 


ecas(^) 





We can rewrite this formula as 


(48) 


ecas(^) = 



2 

— arctan 
7r 




(49) 

Such expressions are known for the total energy, but this local version seems new. 

The local Casimir energy is not, however, equal to the local energy density. 
There are two further terms: 


(i) the energy density contains a term \{d x rf) 2 (where /j(x. t) is the fluctuation 
field), but Casimir energies contain eigenvalues of the field operator which contains 
a term — d 2 ip and our local Casimir energy gets contributions from r)d 2 r ). The 
difference, denoted by Afcas(s), is a double total derivative 


AecasO) = {\d x r]d x ri) - = |< 9 2 <77(37)77(0;)). 


(50) 


s It should also cancel the divergence in the integral / 0 °° x)\ 2 — 1 )7w. This yields 

another amusing formula. 


gj 2^ 2 (s) yjm 2 - u* ^ {x) _ ^ (oo) 

Y j . 2Jm 2 - u?. 5 (<Pk( x ) ~ <Pk (°°))’ 
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= r" I ( ) 


Figure 1: Renormalized tadpoles in the vacuum and in the kink background, respectively. 

The propagator (rj(x)rj(y)) contains a singularity as x tends y, but this singularity 
is ^-independent and cancels due to the space derivatives. So Aecas(^) is a finite 
and smooth function; 

(ii) near the kink, the propagators of // are deformed: they become (complicated) 
expressions for propagation in a kink background. Thus the cancellation of tadpoles 
which we imposed in flat space and which gave us the mass renormalization /ig = 
/r 2 + d/i 2 , no longer holds in the vicinity of the kink. Instead, one has in the kink 
sector 

<p(x,t) = <p K: (x,t) + ifi! (x,t) + rj(x,t), (51) 

where (?/(x, t )) = 0 by definition, and tp\ ( x ) is a mean field induced by the ki nk [Q], 
This mean field gives another contribution to the energy density which we denote 

by Ae( (pi )(x) and which follows from expanding (£(x)) = \ (d x <pd x <p) + \ (U 2 ), 

Ae (v3l) (x) = d x y x d x y K + (i(7 2 )Vi = d x {ip\d x ^p K )- (52) 

The field tpi(x) follows from the vanishing of the expectation value of the field 
equation of the Heisenberg fields $. Using (r)(x, t)) — 0 one easily obtains 

(-5?$ + - (it/ 2 )') = 0 (53) 

= 8%V i - (^t/ 2 )'Vi - i|(if/ 2 )'"(>) 2 ) - i Srn‘ip K (54) 

The sum of the last two terms is again smooth and finite, and if we rewrite this 
equation as 

¥TO) = [dl - ~ l {v 2 (%) - y 2 (oc))3\<p K (x) (55) 

we recognise the Feynman graphs we depicted in Fig. [I]. 

The solution for Lp±(x) is of the form Axrj 0 (x ) + Bd x 7j 0 (x ) (note that (r/ 2 (x)) — 
(rj 2 ( oo)) contains terms with cosfT 4 (mx/2) and rj 0 ~ cosh” 2 (mx/2), and the fluc¬ 
tuation operator d 2 — (\U 2 )" vanishes on r/ 0 ). The term Axr] 0 (x ) can also be written 
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as proportional to <p K (x(l + A)) because rj 0 (x) ~ d x (p K , and this rescaling of x can 
also be written as a rescaling of m (since ipx depends only on -jrix) and a counter 
rescaling of A to keep the prefactor 11 /y/\ invariant: 

d d 

Axrj 0 (x) ~ A ( m Q^ + 2A—^-(rc). (56) 

One can also write this as ^(m, A, x) + Axr] 0 (x ) = ^(m, A, x) where we dis¬ 
cover that this rescaling of the renormalized mass m yields the pole mass rn [§! We 
have not been able to give a similarly simple physical explanation of the rescaled 
coupling A = (m/m) 2 A. 

One can now substitute all expressions to get explicit formulas for the complete 
energy density S(x) for the ki nk (or for any other 1+1 dimensional soliton). One 
can also repeat this exercise for the supersymmetric kink (in this case the only 
difference for ip\ is a different result for m and A, but the term denoted by Bd x r ) 0 is 
the same). However, at this point we refer the reader to the original articles 0.0. 

The local central-charge density has been separately calculated for the susy case 
in [Q] using higher derivative regularization, and also by using susy to transform 
the 7 ■ j anomaly to the sector with the central charge. The explicit result for the 
local central-charge density of [Q] agrees completely with the explicit result for 
the energy density of [0] (where also the explicit local energy density for the non- 
susy case is obtained). In [^9j] a calculation of the integrated central charge can be 
found in global mode regularization, with one cut-off for the Dirac delta function 
in the canonical equal-time (anti)commutation relations and another cut-off for the 
propagators; it is argued that these cutoffs should be the same and this indeed yields 
the correct result. In (§] the anomaly in the local central charge was obtained by 
starting from the definition 

C(y) = j dx8(x - y){\v'{x)U(y) + \tp' (y)U (x)}, (57) 

and not setting x — y too soon. There is a 1 /{x — y) singularity in the propagator 
(■ t]'(x)ri(y )), and expanding x around y in the remaining terms, one finds a finite 
(x — y)/(x — y) term which yields the anomaly. 

So, in conclusion, the nonvanishing one-loop result for the energy density and 
total mass of the minimally susy kink as well as the associated nontrivial modi¬ 
fication of the central charge (density) have been established in the various regu¬ 
larization methods. The specific subtleties of the different methods are now well 
understood, and the origin of the anomalous contribution to the central charge in 
each method clarified, which in particular in dimensional regularization and reduc¬ 
tion shows most interesting facets. 
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Appendix: The DHN program applied to fermions 


The celebrated DHN calculation of the one loop corrections M [V) to the kink mass 
[p due to bosonic fluctuations has been repeated in [15] for the fermionic case, 
using exactly the same steps as DHN did for bosons. We present this calculation 
here for two reasons: (i) to convince skeptics that there is indeed a problem with 
the fermions if one straightforwardly (or better: naively) repeats the same steps, and 
(ii) because there are subtleties with the zero modes which can be clearly illustrated 
by this concrete example. One might anticipate trouble by realizing that for super- 
symmetric boundary conditions (where all non-zero bosonic and fermionic modes 
cancel pairwise) the result for Awould be equal to only the counter term SM 
which is divergent. Some physicists still feel uncomfortable with supersymmetry 
and prefer to stick to older “reliable” methods. The following explicit calculation 
should make it clear that these older methods need updating, for example along the 
lines suggested in the text. 

The field equations for the fermions (= for the fermionic fluctuations) read 
jdip + U'vj = 0, where U 1 = In the representation 7 1 = (q^) and 

7 0 = (° ~ 1 ) the Dirac equation reads 


(d x + U’)^ + = 0 , (~d x + U')il> + d t ^ + 



Iterating and setting ■0 ± = u ± (x)e luJt yields 


<9 2 + u; 2 -m 2 (-tanh 2 — - 
* v 2 2 

r\9 9 9.1 ,9 771X 

dz + or — mr(- tanh-h 

x v 2 2 



0, 

0 . 


(58) 


(59) 

(60) 


The equation for u + is the same as for the bosonic fluctuations 7 (</? = ^k{x) + 
r](x , t)), and hence before imposing boundary conditions, the solutions for u + are 
the same as for //. Given a solution for u + with x / 0, the Dirac equation gives the 
corresponding solution for u~. From the shape of the potentials for the fluctuations 
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(see Fig. JZj), it is clear that 77 and ii) + have a zero mode (a normalizable solution of the 
linearized equations for the fluctuations) but u~ has no normalizable u = 0 solution 
on —00 < x < 00 . However, enclosing the system in a box —L/2 < x < L/ 2, also 
u~ has a zero mode. 

For the bosonic fluctuations the zero mode 770 ~ cosh 2 ( n f-) with strictly ujq = 
0 does not satisfy periodic boundary conditions because its derivative is odd in x, 
but by slightly increasing the energy u, we can achieve that its derivative vanishes 
at the boundaries. Hence, in the bosonic sector there is one almost-zero mode 
with cug O.fi In the second-quantized expression for 77 one finds then a term 
( 2 tu 0 )^ 1 / 2 (a 0 i] 0 (x)e~ luJot + h.c.) which appears on a par with the genuine non-zero 
modes, and hence the almost-zero mode should correspond to one term in the sum 
over zero point energies, just as DHN assumed. 

Imposing even boundary conditions on the fermions^ 

i/j + (~L/ 2) = ^(L/ 2) , r {-L/2) = r(L/ 2), (61) 

we find the following mode solutions for u; 7 ^ 0 : 

= cos (kx ± |<5 + (fc)) cost ot and = — sin (/ex ± ^5~(k)) sine ut, (62) 

7 /) + = — sin (kx ± ^5 + (k)) smut and = cos (kx ± \8~{k )) coscut, (63) 

with ± being a + sign for large positive x and a — sign for large negative x. The 
Dirac equation is satisfied for k > 0 if 

£ = cos \0(h), f = -sm\e{k), 5~(k) =5 + (k) + 9(k), (64) 

9 The other solution for the bosonic fluctuations with uj = 0, given by 
mcosh~ 2 ^ f ' dy cosh 4 4^ does not contribute, even though it is normalizable in the box, 
because it is odd in x and does not tend to zero for large |x|. Hence one cannot make it periodic by 
slightly increasing to 2 . 

1(, Even boundary conditions are not periodic: the derivatives satisfy Robin boundary conditions 

(d x ~ m)ijj + (—L/2) = (d x +m)il> + (L/2) and (d x +m)-i/> _ (— L/2) = (d x — m)ip~(L/ 2) because 
the mass term ?ntanh ^ switches sign between — L/2 and L/2. 
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Figure 3: The phase shift functions 5 + = 5 K , 9, and 5 = 5 + + 9. 




Figure 4: The fermionic zero modes Uq(x) = 1/ cosli 2 (mx/2). u 0 = 0 (concentrated at 
the kink) and u,[ = 0, = cosh 2 (mx/2) (concentrated at the boundary, if any). 


where 5 + {k) = S^ik) is the phase shift of the bosonic fluctuations .[[] The solutions 
with k < 0 are obtained from the solutions with k > 0 by dropping the two minus 
signs in ( |62| ) and t[63j). but for k < 0 ( |54| ) becomes 5~ (k) = 5 + {k ) + 9{k ) + 7r and 
thus the solutions with k < 0 are the same as with k > 0. The cosines satisfy the 
boundary conditions, but the sines must vanish at the boundaries. This yields two 
sets of quantization conditions on k > 0: 


k+L + 5 + (k+) = 2nn + , k n L + 5 (k n ) = n~ 


(65) 


Given the shape of the phase shifts (see Fig. j|) it is clear that n~ = 1,2,3,..., but 
n + = 2,3,4,..., because the solution with n + = 1 (yielding k + = 0) does not 
satisfy the boundary conditions.^ 

We now turn to a closer study of the fermionic zero modes. They are both even, 
one concentrated near the kink and the other near the boundaries (see Fig. |4j). Zero 
modes often correspond to symmetries of the classical action, but the fermionic 
zero mode which is concentrated on the boundary is an example of a zero mode 
which does not correspond to a symmetry, as one might expect since it ceases to be 
a normalizable zero mode on the infinite interval. 

"its explicit form 5pc(k) = —2 arctan(3 mk/ (to 2 — 2k 2 )) is not needed at this point. 

12 The solution with k + = 0 reads u + ~ 3tanh 2 r -kf- — 1, and is even, but u~ = d x + U')u + 
is odd and does not vanish for large |x|. 
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The mode expansions of the fermions read 


- 0 + = 
i>~ = 



( 66 ) 

(67) 


where the sum over n runs over both sets in ( 65]), and where u±(x) are normalized 
to unit-strength plane waves exp i(k n x ± \8 ± (k n )) for large ±|x|. Imposing the 
equal-time canonical anti-commutation relations {^ ± (x,t),^ ± (y,t)} = 8{x — y) 
one finds 

+ ^r{c 0 ,c 0 }u^ (s)< (y) =8{x-y), (68) 


and a similar relation for -0“. To determine the value of the mode anti-commutators, 
we need a completeness relation for the mode functions ?./+ (x), u ( j (x). We go back 
to the second-order differential equation for ip + , and imposing a second boundary 
condition which follows from the Dirac equation 

1> + {~L/ 2) = 2), (d x - m)^ + (-L/ 2) = (d x + m)^ + (L/ 2), (69) 


we obtain a bona-fide selfadjoint elliptic differential operator (with bosonic mode 
operators a n and a*), whose spectrum consists of w+(x), Uq(x). This proves the 
completeness relation 

J2i u n( x ) u n(y) + i u o(xK(y) = s i x ~ y)- ( 7 °) 

n> 0 

Comparing with (|58]) we deduce 

{Cmi ^m,n j {0)> Q)} !• (71) 


The Hamiltonian density for fermions 

H = + U')^ = l -{^ + d t ^ + + (72) 

yields the expected negative contribution to the zero-point energy, (H) = 

~\^Yh n >o u n, and in the density (^ip + (x)i.p + (x)) the zero mode contributes a term 
1jUq(x)uq(x). Due to this factor the two zero modes of the fermions in a box 
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with even boundary conditions contribute one term to the sum over zero-point en¬ 
ergies, just as for the bosonic case, and just as implicitly assumed by DHN.[j] 

We can now compute A. The fermionic modes cancel half of the bosonic 
modes for n > 2, but the bosonic n — 1 mode is left. The bound states and the zero 
modes cancel between bosons and fermions. This yields 


A M (1) 


_ l52 U ( k n ) + SM 


n>l 




n> 2 


I Yp^-^+SM 


n> 2 
poo 


\ m + \ 


dk 

dk dcu . r , , 7T — 2 
2^*1 + 


(73) 


where we used that = k~ + 6{k~)/L. This is the correct answer to an incorrect 
question, because this value for M ^ contains spurious boundary energy. In the text 
it is discussed how to separate off the spurious boundary energy, and the correct 
result is 

M^ = ~—. (74) 

Ztt 

If one repeats the same calculations for the sine-Gordon system, one can compare 
with the exact result obtained from the Yang-Baxter equation, and finds indeed that 
these two results agree after removing the boundary energy. 
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